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I. 


INTRODUCTION 


A.  ANALYSIS  OBJECTIVES 

This  thesis  implements  a  model  for  the  transmission  of  acoustic  energy  in  shallow 
maritime  waters.  High-frequency  (40-70  kHz)  sound  is  propagated  numerically  through 
shallow  (5-40  m)  water  channels.  Parametric  analysis  determines  the  influence  of 
environmental  conditions  by  varying  a  specific  factor  while  holding  all  other  factors 
constant.  This  parametric  analysis  includes  changes  in  the  sound-speed  profile  (SSP) 
gradient,  changes  in  transmission  frequency,  changes  in  transmission  range,  changes  in 
channel  depth,  changes  in  the  bottom  sediment  material  of  the  channel,  and  changes  in 
the  height  of  the  transmitter  and  the  receiver  in  the  channel. 

Finally,  the  lessons  learned  from  the  parametric  analyses  are  used  to  characterize 
two  maritime  environments  where  experimentation  is  being  conducted  in  2007.  The 
testing  sites  are  St.  Andrew  Bay,  in  Panama  City,  FL,  and  Chesapeake  Bay,  in  Little 
Creek,  VA.  The  purpose  of  characterizing  these  environments  is  to  predict  the 
performance  of  experimental  communication  systems  before  they  are  deployed. 

B.  UNDERWATER  COMMUNICATIONS  AND  NETWORKS 

Seaweb  is  an  experimental  undersea  network  technology  consisting  of  medium 
range  (500-5000  m)  links  between  network  nodes.  The  Seaweb  wide-area  network 
(WAN)  concept  includes  local-area  networks  (LANs)  having  a  sophisticated  central  node 
that  collects  and  integrates  undersea  data  from  a  set  of  simpler  peripheral  nodes.  This 
local-area  network  with  centralized  topology  is  called  Seastar  and  is  motivated  by  the 
desire  for  an  additional  tier  of  local-area  communications  compatible  with  and 
complementary  to  Seaweb  wide-area  acoustic  communications.  Seastar  uses  a  higher- 
frequency  portion  of  the  acoustic  spectrum  matched  to  the  short  range  (50-500m)  links. 
The  high  frequencies  allow  for  a  larger  spectral  bandwidth  and  thus  a  higher  rate  of 
information  transfer  in  the  LAN.  The  Seastar  LAN  and  its  place  in  the  Seaweb  WAN  are 
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illustrated  in  Figure  1,  with  the  red  nodes  representing  peripheral  nodes  and  the  green 
nodes  representing  central  nodes. 


Figure  1.  Seastar  local-area  network  (LAN)  concept. 

The  Seastar  central  node  integrates  data  received  from  the  peripheral  nodes,  and 
in  turn  transmits  fused  information  to  other  central  nodes  or  to  a  gateway  node  for 
transmission  to  command  centers.  Communications  amongst  central  nodes  and  gateway 
nodes  occur  via  Seaweb  networking  in  a  lower  band  of  the  acoustic  spectrum. 

The  importance  of  Seastar  is  that  it  allows  an  order-of-magnitude  increase  in  node 
density  in  the  Seaweb  infrastructure,  and  the  use  of  a  higher-frequency  band  means  that 
central  nodes  can  maintain  both  Seaweb  and  Seastar  communications  simultaneously 
without  mutual  interference.  Seastar  allows  for  high-bandwidth,  low-power 
communications,  advancing  environmental  sensing,  intelligence  preparation  of  the 
battlespace,  maritime  special  ops,  and  littoral  surveillance. 
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II.  PHYSICAL  BASIS  FOR  THE  NUMERICAL  MODEL 


A.  THEORETICAL  BACKGROUND 


The  wave  equation  is  a  partial  differential  equation  describing  the  motion  of  a 
wave  in  a  medium.  The  acoustic  wave  equation  is  a  linear  approximation  retaining  only 
first-order  terms  of  the  wave  equation.  The  propagation  of  sound  in  an  elastic  medium 
(such  as  seawater)  is  described  mathematically  by  solving  the  acoustic  wave  equation 
while  applying  the  appropriate  boundary  and  medium  conditions.  The  acoustic  wave 
equation  relates  acoustic  pressure  p  to  spatial  coordinates  x,  y,  z  and  time  t,  and  may  be 
written  in  the  standard  form  of  the  wave  equation  as 


vV-Tiz= o 

c  dt 


(2.1) 


d  d  d 

where  V  is  the  gradient  operator,  V  =  —  h - h  — ,  and  c  is  the  speed  of  sound  in  the 

dx  dy  dz 

medium  and  may  vary  with  the  spatial  coordinates.1 


1.  Wave  Theory 

Wave  theory  describes  sound  propagation  in  terms  of  characteristic  functions 
called  normal  modes,  each  of  which  is  a  solution  to  Equation  (2.1).  These  normal  modes 
are  combined  additively  to  satisfy  the  given  boundary  and  medium  conditions.  Although 
this  approach  provides  a  formally  complete  solution  for  all  frequencies,  it  has  the 
disadvantages  of  not  presenting  an  intuitive  picture  of  the  distribution  of  radiated  energy, 
not  easily  handling  real  boundary  conditions,  and  requiring  significant  computation  since 
it  tries  to  model  the  entire  wave.2 


1  Robert  Urick,  Principles  of  Underwater  Sound  3rd  Edition  (Los  Altos:  Peninsula  Publishing,  1 983), 
120-122. 

2  Ibid,  122. 
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2.  Ray  Theory 


Ray  theory  is  an  alternative  theoretical  approach  to  solving  the  wave  equation. 
The  main  tenets  of  ray  theory  are  the  existence  of  a  wavefront  along  which  the  phase  of 
the  solution  is  constant,  and  the  existence  of  rays  that  describe  the  spatial  location  of 
acoustical  energy  radiating  from  the  source  in  a  manner  analogous  to  optical  ray  theory. 
Ray  tracing  programs  are  commonly  used  in  underwater  acoustics  for  modeling  the 
propagation  of  high-frequency  acoustic  waves  as  a  function  of  time.3 

Ray  tracing  traditionally  involves  integrating  a  set  of  differential  equations  called 
the  ray  equations,  which  describe  the  ray’s  trajectory.  These  equations  are  governed  by 
given  initial  conditions  in  order  to  trace  the  path  of  a  ray  as  it  propagates  away  from  the 
source.  The  amplitude  of  a  ray  is  determined  by  the  cross-section  of  the  ray  tube  bounded 
by  adjacent  rays. 

The  main  drawback  associated  with  ray  tracing  is  the  existence  of  shadow  zones, 
which  are  zones  through  which  no  rays  pass,  resulting  in  a  pressure  field  of  zero 
everywhere  within  them.  In  actuality,  there  is  always  some  diffraction  of  sound  into  the 
areas  given  as  shadow  zones  in  ray  tracing,  resulting  in  a  discrepancy  between  the  exact 
solution  and  that  predicted  by  ray  tracing  methods.4 

Ray  tracing  is  a  useful  tool  for  modeling  the  propagation  of  high-frequency 
sound.  Ray  tracing  is  intuitive  and  radiated  sound  is  readily  visualized  and  calculated. 

B.  MATHEMATICAL  DERIVATION 

Ray  theory  is  mathematically  derived  from  the  acoustic  wave  equation.  The 
derivation  yields  the  time  delay  of  the  ray  from  source  to  receiver,  the  ratio  of  the 
pressure  amplitude  at  the  desired  location  to  the  pressure  amplitude  at  a  reference  point, 
and  the  transmission  loss.  These  functions  are  used  to  construct  ray  and  impulse  response 
diagrams. 


3  Urick,  122. 

4  Ibid. 
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The  basis  for  ray  theory  is  the  acoustic  wave  equation,  given  in  Equation  (2.1). 
Separation  of  variables  shows  pressure  depends  on  the  three-dimensional  position  vector 
x  =  (x,  y,  z)  and  time  t, 

p  =  P(x)-T(t)  ,  (2.2) 

so  that  we  are  only  separating  the  spatial  and  time  dependences.  Substituting  this  form 
into  Equation  (2.1)  and  taking  the  separation  constant  to  be  k2  we  obtain5 

rI2T 

V2P  +  k2P  =  0,  — T  +  k2c2T  =  0.  (2.3) 

dt 2 

The  first  equation  in  (2.3)  is  the  time-independent  version  of  the  acoustic  wave 
equation,  called  the  Helmholtz  equation.  If  we  substitute  k  =  col  c ,  then 

V2p  +  -^-p  =  0,  (2.4) 

c  (x) 


where  c(x)  is  the  tri-dimensional  sound  speed  and  co  is  the  angular  frequency  of  the 
source.  Jensen  et  al.  develop  a  solution  of  the  Helmholtz  equation  of  the  form6 


/><xW'wX 

7=0 


Mx) 
(icoy  ’ 


(2.5) 


where  r(x)  is  the  time  it  takes  for  the  sound  to  reach  location  x  and  A(x)  is  the  amplitude 
of  the  signal  at  x.  Equation  (2.5)  is  called  the  ray  series.  The  first  and  second  derivatives 
of  the  ray  series,  retaining  only  the  first-order  terms  (a  high-frequency  approximation) 
gives  the  following  infinite  sequence  of  equations  for  the  functions  r(x)  and  A(x), 

0(a>2):  |Vr|2=— T_  (2.6) 

c  (x) 


0(®):  2Vr-V4o+(V2r)4,  =0 


(2.7) 


5  Kenneth  F.  Riley  et  al.,  Mathematical  Methods  for  Physics  and  Engineering  (Cambridge:  Cambridge 
University  Press,  1998),  566. 

6  Finn  B.  Jensen  et  al.,  Computational  Ocean  Acoustics  (New  York:  Springer- Verlag,  2000),  143-151. 
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so  that  the  solution  to  the  Helmholtz  equation  is  a  product  of  a  phase  function,  called  the 
eikonal  equation,  and  an  amplitude  function,  called  the  transport  equation. 


To  solve  the  eikonal  equation,  a  family  of  rays  perpendicular  to  the  wavefronts  of 


x(x)  are  introduced.  This  family  of  rays 
coordinates  in  terms  of  the  arc  length,  s. 
differential  equation 

dx  _ 
ds 


defines  a  new  coordinate  system  in  ray 
The  ray  trajectory  is  then  defined  by  the 

cVr  (2.8) 


which  can  be  squared  to  yield 


dx 


|2 


ds 


(2.9) 


Substituting  Equation  (2.6)  into  Equation  (2.9)  shows  that  the  tangent  vector 
dx/ds  has  unit  length.  Taking  the  Cartesian  coordinates  of  Equation  (2.8)  and 
differentiating  with  respect  to  the  arc  length  s  gives  the  resultant  vector  equation  for  the 
ray  trajectories 


d  f  1  dx 
ds  vc  ds 


(2.10) 


In  cylindrical  coordinates  (r,  z)  the  ray  equations  can  be  written  in  the  first-order 

form 


dr 

ds 


c£(s). 


dz 

ds 


=  cg(s), 


d%  _  1  dc 

ds  c 2  dr 
dg  _  l  dc 
ds  c2  dz  ’ 


(2.11) 


where  integrating  these  equations  gives  the  trajectory  of  the  ray  [r(s),  z(s)].  For  these 
equations,  c  and  C  are  auxiliary  variables  arbitrarily  introduced  so  the  equations  may  be 
written  in  first-order  form. 


The  initial  conditions  for  solving  the  ray  equations  are  that  the  ray  starts  at  the 
source  position  (rs,  zs)  with  a  specified  vertical  launch  angle  6 ,  as  shown  in  Figure  2. 
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Schematic  of2-D  ray  geometry 


Figure  2.  Schematic  of  2-D  ray  geometry. 


Thus  we  have 


r  = 


cos# 

40? 


(2.12) 


z  = 


sin# 

4?' 


(2.13) 


To  obtain  the  pressure  field  the  phase  and  amplitude  of  each  ray  are  required.  The 
phase  is  obtained  by  solving  the  eikonal  equation  in  the  ray  coordinate  system. 
Expanding  Equation  (2.6), 


„  „  1 
Vt  -Vt  =  —  , 
c 

so  that  substituting  the  value  of  V  t  from  Equation  (2.8), 

^  1  dx  1 

Vr--— =  — , 
c  as  c 

then 

dr  _  1 
ds  c 


(2.14) 


(2.15) 


(2.16) 
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Integrating  this  equation  with  respect  to  s  yields 


5  5  1 

\dr  =  - ds 

l  l  c(s) 

t(s)-t(0  )=[-}— ds 
0  c(s) 

r(s)  =  r(0)  +  j-}—ds. 

o  c(s ) 


(2.17) 


The  integral  term  in  this  equation  is  the  travel  time  along  a  ray  (time  delay),  so 
that  the  phase  of  the  wave  is  delayed  in  accordance  with  its  travel  time. 

The  amplitude  is  obtained  by  solving  the  transport  equation.  Since  the  rays  are 
perpendicular  to  the  wavefronts,  Equation  (2.7)  can  be  rewritten  as 

2  dx 


c  ds 


■VA0+(Vzt)A0=  0 


(2.18) 


which  gives 


2  dA0 
c  ds 


+  (V2t)A0  =  0. 


(2.19) 


This  equation  thus  states  that  the  amplitude  along  a  ray  changes  in  relation  to  the 
spreading  of  a  ray  tube,  which  is  the  space  bounded  by  two  adjacent  rays,  the  solution  to 
the  transport  equation  can  be  simplified  by  using  the  three-dimensional  Jacobian 
determinant .  Since  the  Jacobian  satisfies  the  following  condition, 

O' 

J  ds  I 


V2r  =  — — 


(2.20) 


\cj 


Equation  2.19  then  becomes 

_  dAn 


ds 


-  + 


c  d 
J  ds 


(-) 

UJ 


Ao=0. 


(2.21) 


Integrating  this  equation  from  0  to  s  gives  the  solution  for  the  transport  equation 
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A0(s)  =  A0(  0) 


c(s)J(0) 

c(0)J(s) 


(2.22) 


The  sound  pressure  for  a  point  source  in  a  homogenous  medium  is  defined  as 


p0(s)  =  A0(s)e^  = 


(2.23) 


So  it  follows  that  Aq(s)  and  to(s)  are  given  by 


(2.24) 


r0O)  = 


(2.25) 


Substituting  the  limit  of  the  quantity  ^o(O)  •  /(0)  into  Equation  (2.22)  yields 


Ms)  =  — 


c(s ) cos  6 


An  c(0)/(s) 


(2.26) 


so  the  sound  pressure  as  a  function  of  arc  length  s  is 


,  ,  1  c(5)cos^ 

pis)  = - e  0 

An  c(0)J(s) 


(2.27) 


Thus,  the  pressure  field  is  obtained  by  dividing  the  energy  of  the  point  source 
among  each  of  the  ray  tubes.  Since  a  point  source  in  a  shallow  water  channel  follows 
spherical  spreading,  the  sound  pressure  amplitude  of  a  ray  tube  diminishes  with  range  5 
by  the  factor  Ms  .  The  ratio  of  the  pressure  at  a  point  a  distance  s  from  the  source,/?^),  to 
the  intensity  measured  at  1  m  from  the  source,  po,  defines  the  transmission  loss  at  that 
point.  For  convenience,  the  transmission  loss  along  each  ray  tube  is  usually  given  in 
decibels  (dB): 


TL(s)  =  -20  log - 


(2.28) 
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The  total  pressure  at  the  receiver  will  depend  on  the  way  in  which  the  rays 
constructively  and  destructively  interfere  with  each  other. 
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III.  BELLHOP 


A.  MODEL  DESCRIPTION 

Bellhop  is  a  Gaussian  beam  tracing  acoustic  propagation  model  developed  in 
1987  by  Porter  and  Bucker7  at  the  Space  and  Naval  Warfare  Systems  Center  in  San 
Diego.  Bellhop  was  chosen  for  this  analysis  since  it  has  proven  to  be  an  accurate 
modeling  tool  for  high-frequency  (>1  kHz)  transmissions. 

Bellhop  is  a  range-independent  program,  which  means  that  it  assumes  that  the 
sound-speed  profile  (SSP),  the  surface  and  bottom  properties,  and  the  bathymetry  of  the 
ocean  floor  do  not  vary  with  range.  Propagation  is  computed  for  a  2-dimensional  vertical 
slice  of  the  medium,  under  the  assumption  that  this  result  applies  for  all  bearings.  These 
assumptions  allow  a  reduction  in  the  computational  times.  Bellhop  includes  an  option 
that  allows  the  program  to  import  a  bathymetry  file,  but  this  function  was  not  used  as  part 
of  the  following  set  of  analysis. 

Bellhop  uses  Gaussian  beam  tracing  to  produce  the  ray-trace  and  impulse- 
response  plots.  Gaussian  beam  tracing  uses  standard  ray  tracing  to  describe  a  central  ray, 
about  which  a  beam  is  formed  whose  intensity  decays  on  a  normal  from  the  ray  following 
a  Gaussian  distribution.  Beam  tracing  is  computationally  more  efficient  while  being 
mathematically  equivalent  to  standard  ray  tracing  and  it  allows  for  a  transition  of  the 
beams  into  shadow  zones,  thus  mitigating  an  aforementioned  drawback  of  simple  ray 
tracing. 

Bellhop  is  able  to  trace  a  user-defined  fan  of  rays  emitted  from  a  sound  source 
whose  coordinates  and  frequency  are  specified  by  the  user.  It  also  calculates  the 
amplitude  and  arrival  time  of  an  impulse  signal  arriving  at  the  receiver  via  each  ray  path. 
Bellhop  can  also  compute  the  intensity  of  the  sound  field  at  the  receiver  by  finding  only 
the  rays  that  reach  the  receiver.  These  rays  are  called  eigenrays. 

7  Michael  B.  Porter  and  Homer  P.  Bucker,  “Gaussian  beam  tracing  for  computing  ocean  acoustic 
fields,”  The  Journal  of  the  Acoustical  Society  of  America  82  (1987):  1349-1359. 
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B.  INPUT  ENVIRONMENT  FILE 

The  main  input  file  to  Bellhop  is  called  the  environment  file.  The  environment  file 
is  a  simple  text  file  given  in  a  structured  format  specifying  the  details  of  the  environment, 
the  source  and  receiver  characteristics,  and  the  type  of  analysis  to  be  performed.  The 
environment  information  includes  the  sound-speed  profile,  the  depth  of  the  channel,  the 
channel  bottom  composition,  and  surface  boundary  conditions.  The  source  and  receiver 
characteristics  include  the  transmission  frequency,  the  source  and  receiver  locations,  the 
angle  limits  for  the  fan  of  beams  and  the  number  of  beams  used.  The  type  of  analysis 
specified  dictates  the  output  of  the  algorithm  in  the  form  of  a  ray  file,  an  eigenray  file,  an 
amplitude-delay  file,  or  any  of  three  different  types  of  transmission  loss  files  (coherent, 
incoherent  and  semi-coherent). 

C.  OUTPUT  PRODUCTS 

1.  Ray-Trace  Plots 

Information  in  the  output  ray  file  can  produce  a  plot  showing  the  paths  traveled  by 
all  the  rays  specified  in  the  environment  file.  The  plot  has  depth  on  the  vertical  axis  and 
range  on  the  horizontal  axis.  The  rays  in  the  plots  presented  in  Appendix  B  are  in  three 
colors.  A  red  ray  is  a  direct  path  between  the  source  and  the  receiver,  generally  signifying 
that  this  is  the  strongest  path.  A  blue  ray  has  only  one  reflection  between  the  source  and 
receiver,  which  signifies  an  intermediate  strength  path.  A  black  ray  has  multiple  surface 
and  bottom  reflections  between  source  and  receiver,  meaning  this  is  the  weakest  type  of 
path.  Ray-trace  plots  show  all  rays  regardless  of  whether  they  pass  through  the  receiver 
or  not.  This  is  useful  for  visualizing  the  radiated  energy  distribution  in  the  channel. 

2.  Eigenray  Plots 

Bellhop  can  selectively  plot  the  paths  traveled  by  only  those  rays  that  reach  the 
receiver.  Since  Bellhop  uses  geometric  beams  in  its  calculations,  the  eigenrays  include 
those  ray  tubes  where  the  receiver  is  between  the  adjacent  rays.  This  plot  has  the  same 
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axes  and  color  scheme  as  the  ray-trace  plot.  This  plot  displays  the  multipath  propagation 
from  source  to  receiver. 

3.  Impulse  Response  Plots 

Bellhop  can  also  create  an  amplitude-delay  file,  which  is  a  text  file  stating 
transmission  frequency,  channel  depth,  source  and  receiver  geometry  and  number  of 
eigenrays.  Also  included  is  information  about  each  eigenray  such  as  amplitude,  time 
delay,  phase  shift  at  the  receiver,  and  the  number  of  surface  and  bottom  bounces 
experienced  by  the  eigenray.  An  algorithm  using  the  information  in  this  file  to  plot  the 
impulse  response  h(z)  is  developed  and  listed  in  Appendix  A.  The  impulse  response  h(z) 
is  defined  as  the  output  of  the  receiver  to  a  very  brief  signal  from  the  source,  called  an 
impulse.8  The  plot  has  the  dimensionless  amplitude  at  the  receiver  in  the  vertical  axis  and 
the  time  delay  in  seconds  in  the  horizontal  axis.  The  impulse  response  plot  shows  the 
time  dispersion  of  an  initial  impulse  at  the  source  into  multiple  impulses  at  the  receiver. 
This  time  dispersion  is  caused  by  the  different  pathlengths  and  travel  times  associated 
with  each  ray.  An  important  implication  of  this  multipath  time  dispersion  of  signal  energy 
is  the  possibility  of  intersymbol  interference  for  acoustic  communications. 

The  implementation  of  Bellhop  used  here  neglects  the  possibility  of  dispersion  in 
the  frequency  domain,  such  as  that  occurring  as  a  result  of  Doppler  shift  and  Doppler 
spread  associated  with  a  moving  source  and/or  receiver,  and  a  non- stationary  channel. 


8  Johnny  R.  Johnson,  Introduction  to  Digital  Signal  Processing  (Englewood  Cliffs:  Prentice  Hall, 
1989),  51. 
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IV.  SHALLOW  WATER  PARAMETRIC  ANALYSIS 


A.  STANDARD  CHANNEL  DEFINITION 

The  standard  channel  is  defined  for  this  parametric  analysis  to  be  10  m  deep  with 
a  sand  bottom  boundary  and  a  vacuum  over  the  surface  boundary,  a  source-to-receiver 
horizontal  range  of  500  m,  and  both  source  and  receiver  1  m  above  the  seabed.  The 
standard  transmission  is  assumed  to  be  at  70  kHz  with  source  power  of  1  W.  Three 
generic  sound-speed  profiles  are  analyzed.  These  three  profiles  (salinity-driven, 
temperature-driven  and  isospeed)  represent  typical  conditions  encountered  in  shallow 
water  and  are  illustrated  in  Figure  3. 

Salinity-Driven  Environment  Isospeed  SSP  Temperature-Driven  Environment 


Figure  3.  Sound-speed  profiles  (SSP)  used  in  the  parametric  analyses 
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B.  TRANSMITTER 


The  source  level  ( SL )  of  a  transmitter  is  defined  as  9 

££  _  source @\m  (4  1) 

I  ref 

The  source  level  of  an  omnidirectional  projector  is  referenced  to  a  standard  range 
of  1  m  from  its  acoustic  center.  At  a  radius  of  1  m,  the  acoustic  center  is  surrounded  by  a 
sphere  of  surface  area 

APh  =  4^2  =>  Asph  =  12.6  m2.  (4.2) 

If  the  power  of  the  transmitter  is  given  by  Pxmt,  the  source  intensity  at  1  m  is  given 
by 

/,  =  — W/m2.  (4.3) 

12.6 

The  standard  reference  pressure  used  in  underwater  acoustics  is  1  pPa,  which  is 
approximately  equivalent  to  0.67-10  18  W/m2.10  Using  this  conversion  and  Equation 
(4.3)  the  source  level  becomes 

SL  =  KmL - I —  =>  SL  =  10 log (Pxmt )  + 1 70.8  dB  re  1  pPa  @  1  m.  (4.4) 

12.6  0.67 -10'18  xmt 

If  we  assume  the  source  has  a  power  of  1  W,  then  the  source  level  is 

SL  =  101og(l)  +  170.8dB  re  1  pPa  @  1  m  =>  SL  =  170.8  dB  re  1  pPa  @  1  m.  (4.5) 

The  amplitude-delay  file  produced  by  Bellhop  contains  the  ratio  of  the  pressure 
amplitude  at  the  desired  distance  divided  by  the  pressure  amplitude  at  the  source.  The 
transmission  loss  is  then  calculated  by  using  Equation  2.28.  The  amplitude  of  the  sound 
signal  at  the  receiver  is  obtained  by 


9  Ashley  D.  Waite,  SONAR  for  Practising  Engineers  (West  Sussex:  John  Wiley  &  Sons  Ltd.,  2002),  4- 
5.  '  ' . . 

Lawrence  E.  Kinsler  et  al.,  Fundamentals  of  Acoustics  4th  Edition  (New  York:  John  Wiley  &  Sons, 
2000),  131. 
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A(s)  =  SL-TL(s) , 


(4.7) 


where  A(s j  is  given  in  dB.  This  amplitude  is  used  as  the  vertical  axes  for  all  the  summary 
figures  and  impulse  response  figures  developed  in  this  analysis,  since  this  amplitude  is  a 
function  of  all  the  parameters  being  varied. 


C.  MARSH  &  SCHULKIN  EQUATION 


In  1962,  Marsh  and  Schulkin11  developed  an  equation  for  transmission  loss  in 
shallow  water.  The  transmission  loss  in  these  media  is  dependent  upon  variables  of  the 
sea  surface,  water  channel  and  sea  botton.  Since  the  channel  is  very  sensitive  to  changes 
of  these  variables,  the  predicted  transmission  loss  can  only  be  approximated.  The 
empirical  basis  for  the  Marsh  &  Schulkin  equation  is  100,000  sound  level  measurements 
they  made  in  shallow  water.  The  applicability  of  the  equation  is  contingent  upon 
determination  of  a  factor  H,  which  is  defined  as 


H  = 


^3  3  A 

D 


1/2 


8 


j 


(4.8) 


where  D  is  the  water  depth  in  meters.  Table  1  shows  the  value  of  H  for  the  different 
channel  depths  used  in  this  study. 


D 

H 

(m) 

(km) 

5 

1.44 

10 

2.03 

20 

2.87 

40 

4.06 

Table  1 .  Marsh  &  Schulkin  equation  parameter  H  for  different  channel  depths 


There  are  three  different  transmission  loss  equations  depending  upon  the 
relationship  between  H  and  the  transmission  range  r.  For  this  analysis  only  one  equation 


H.  W.  Marsh  et  al.,  “Shallow  Water  Transmission.”  The  Journal  of  the  Acoustical  Society  of 
America  34  (1962):  863-864. 
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will  be  used  since  all  of  the  transmission  ranges  used  in  this  analysis  are  smaller  than  H. 
For  this  condition,  the  transmission  loss  in  dB  is  given  by 


71  =  20*  log 


r 

dm 


+  a 


r 

V14 


+  60  dB, 


(4.9) 


where  r  is  the  transmission  range  in  km,  a  is  the  absorption  coefficient,  1/0.914  is  the 
conversion  factor  between  km  and  kyd,  and  60  dB  accounts  for  the  factor  of  1000 
between  m  and  km.  Equation  (4.9)  is  different  from  Equation  (2.28)  since  the  former 
computes  transmission  loss  on  the  channel  while  the  latter  computes  the  transmission  loss 
on  a  single  ray  tube.  The  Marsh  &  Schulkin  equation  for  transmission  loss  is  used  as  a 
baseline  for  comparison  with  the  transmissions  in  this  analysis. 


D.  SOUND-SPEED  PROFILE  (SSP) 

Variations  in  sound  velocity  result  in  refraction  of  the  sound  rays.  Sound  velocity 
is  influenced  by  temperature,  salinity  and  water  pressure.  Temperature  effects  are  mainly 
due  to  surface  solar  heating,  nocturnal  cooling  and  wind  action.  The  salinity  of  seawater 
can  vary  due  to  precipitation,  fresh-water  runoff  and  mixing.  Water  pressure  is  a  linear 
function  of  depth.  The  combination  of  these  factors  can  lead  to  a  wide  array  of  sound- 
speed  profiles.  The  purpose  of  this  section  is  to  analyze  seven  different  profiles  to 
identify  advantageous  and  disadvantageous  conditions  for  transmission.  The  salinity- 
driven  and  temperature-driven  profiles  are  plotted  in  Figure  3,  and  the  five  linear  sound- 
speed  profiles  are  plotted  in  Figure  4. 


-2  SVP  Gradient 


+1  SVP  Gradient 


Figure  4.  Linear  sound-speed  profiles  (SSP) 
18 


Figure  5  shows  the  amplitudes  obtained  for  seven  different  types  of  sound 
velocity  profiles.  The  isospeed  and  temperature-driven  profiles  are  observed  to  be  similar 
in  strength  since  they  both  have  direct  paths  between  the  transmitter  and  receiver.  The  -2 
SSP  gradient  profile  is  observed  to  strongly  focus  the  sound  energy  towards  the  bottom  of 
the  channel.  It  is  hypothesized  that  this  focusing  of  energy  and  constructive  interference 
contribute  to  this  profile  being  the  strongest;  conversely,  it  is  hypothesized  that  the  +2 
SSP  gradient,  a  strongly  linear  upward  refracting  profile  focuses  the  energy  towards  the 
top  of  the  channel  and  therefore  away  from  the  receiver,  so  it  displays  the  smallest 
amplitude.  Small  linear  changes  in  the  profile  (+1  and  -1  SSP  gradients)  and  the  salinity- 
driven  profile  are  observed  to  disrupt  the  direct  path  of  transmission  but  not  the 
secondary  paths,  so  they  reduce  the  magnitude  of  the  largest  amplitude  but  not  as 
drastically  as  the  +2  SSP  gradient  profile. 


Sound-Velocity  Profile 

Figure  5.  Impact  of  sound-speed  profile  (SSP) 
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E.  TRANSMISSION  FREQUENCY 


When  sound  propagates  in  the  ocean,  part  of  the  acoustic  energy  is  absorbed  and 
scattered  by  the  medium.  The  combination  of  absorption  and  scattering  effects  is  called 
attenuation  and  it  is  frequency-dependant.  The  quantity  a  is  the  logarithmic  absorption 
coefficient,  and  it  is  expressed  in  units  of  dB/km.  The  following  expression,  developed  by 
W.  H.  Thorp,  relates  frequency  in  kHz  to  the  attenuation  coefficient12 

oil  f2  44  f 2 

a  ^  3.3x1 0“3  +  '  J-  + - J- — 7  + 3. Oxl 0”4/2 .  (4.10) 

1  +  /2  4100  +  f2 

The  attenuation  coefficient  a  for  frequencies  in  the  40-70  kHz  band  is  calculated 
using  Equation  (4.10)  and  presented  in  Table  2  and  Figure  6.  These  results  show  that  as 
the  frequency  increases,  the  value  of  the  attenuation  coefficient  increases  as  well,  which 
leads  to  a  greater  transmission  loss  and  thus  smaller  amplitude  at  the  receiver.  The  theory 
thus  predicts  a  reduction  in  amplitude  at  the  receiver  for  each  increase  in  frequency. 


Frequency  (kHz) 

a  (dB/km) 

40 

12.5 

45 

14.7 

50 

16.8 

55 

18.8 

60 

20.7 

65 

22.5 

70 

24.1 

Table  2.  Attenuation  coefficients  for  the  40-70  kHz  band 


12  Urick,  108. 
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Figure  6.  Attenuation  coefficient  vs.  frequency 

The  results  for  the  frequency  parametric  analysis  are  summarized  in  Figure  7.  The 
figure  shows  that  the  isospeed  environment  and  the  temperature-driven  environment 
behave  as  predicted  by  the  theory,  with  diminishing  amplitudes  at  the  receiver  for 
increasing  frequencies.  The  figure  also  shows  that  the  profiles  behave  similarly  to  the 
Marsh  &  Schulkin  equation  prediction. 

In  the  salinity-driven  case,  the  40-45  kHz  band  and  the  50-70  kHz  band  also 
behave  as  predicted.  In  the  45-50  kHz  band,  however,  there  is  an  increase  in  the 
amplitude  at  the  receiver  for  increasing  frequency.  The  reason  for  this  anomaly  is  that  the 
Bellhop  algorithm  apparently  attributes  the  maximum  impulse  for  the  40-50  kHz  band 
and  the  50-70  kHz  band  to  different  paths.  Since  the  frequency  is  the  only  parameter 
being  changed,  the  geometrical  distribution  of  the  rays  should  not  change  and  the  path  or 
the  maximum  impulse  should  be  the  same  for  all  frequencies. 
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Figure  7.  Impact  of  transmission  frequency 


F.  TRANSMISSION  RANGE 

Sound  produced  by  an  omnidirectional  point  source  in  a  homogenous,  lossless 
and  unbounded  medium  is  radiated  equally  in  all  directions  so  as  to  be  equally  distributed 
over  the  surface  of  a  sphere  surrounding  the  source.13  The  intensity  decreases  as  the 
square  of  the  range,  so  the  transmission  loss  is  then 

TL  =  20  log(r)  +  a-(r/l  000) ,  (4.11) 

where  r  is  in  meters,  a  is  the  attenuation  coefficient  given  in  dB/km  and  the  transmission 
loss  has  units  of  dB. 

Sound  produced  by  a  point  source  propagating  through  a  medium  limited  by 
plane-parallel  surface  and  bottom  boundaries  is  approximated  as  distributed  over  the 
surface  of  a  cylinder  having  a  radius  equal  to  range  and  a  height  equal  to  the  depth  of  the 


13  Urick,  100 
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channel.14  The  intensity  of  the  sound  is  then  inversely  proportional  to  the  surface  of  the 
cylinder,  and  the  transmission  loss  increases  proportionally  to  the  transmission  range. 

TL  =  101og(r) +  «•(>/ 1000) ,  (4.12) 

where  r  is  in  meters,  a  is  the  attenuation  coefficient  given  in  dB/km,  and  the  transmission 
loss  has  units  of  dB. 

Figure  8  is  a  comparison  of  the  amplitude  at  the  receiver  for  a  1-W,  70-kHz 
source  for  three  different  transmission  loss  assumptions  (cylindrical,  spherical,  and  Marsh 
&  Schulkin  equation).  The  figure  shows  that  cylindrical  spreading  causes  the  least 
amount  of  transmission  loss,  since  the  energy  is  reflected  back  into  the  channel  from  the 
boundaries.  Spherical  spreading  causes  more  severe  transmission  loss  since  no  energy  is 
reflected  back  into  the  channel.  The  Marsh  &  Schulkin  equation  produces  the  greatest 
amount  of  transmission  loss  and  most  closely  matches  the  Bellhop  results.  All  three 
transmission  loss  equations  display  similar  trends  for  increasing  range. 


Figure  8.  Amplitude  vs.  range  for  three  different  types  of  transmission  loss 


14  Ibid,  102. 
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Figure  9  shows  the  amplitude  at  the  receiver  for  the  three  different  types  of  sound- 
speed  profile  at  different  transmission  ranges.  The  figure  shows  that  the  temperature- 
driven  profile  gives  the  largest  amplitude  for  distances  less  than  475  m.  For  distances 
greater  than  475  m,  the  isospeed  profile  gives  the  largest  amplitude.  The  figure  also 
shows  that  the  temperature-driven  and  salinity-driven  environments  follow  similar  trends, 
and  that  the  isospeed  profile  closely  resembles  the  prediction  from  the  Marsh  &  Schulkin 
equation. 


Figure  9.  Impact  of  transmission  range 


G.  CHANNEL  DEPTH 

Four  different  channel  depths  (5,  10,  20  and  40  m)  are  analyzed.  Table  1  shows 
that  these  depths  fall  within  the  condition  of  the  Marsh  &  Schulkin  equation.  To  obtain 
the  5  m  deep  channel,  the  10  m  standard  profiles  are  cut  at  the  5  m  mark.  To  obtain  the  20 
m  and  40  m  deep  channels,  it  is  assumed  that  all  the  profiles  are  nearly  isospeed  at  the 
bottom.  Thus,  isospeed  extrapolation  completes  the  profile.  All  tests  have  the  transmitter 
and  receiver  1  m  from  the  bottom  of  the  channel.  Figure  10  shows  the  amplitude  at  the 
receiver  for  the  three  sound-speed  profiles  at  different  channel  depths.  The  figure  shows 
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that,  since  there  is  no  variation  in  the  Marsh  &  Schulkin  equation  parameters,  the 
amplitude  at  the  receiver  is  predicted  to  be  constant.  The  isospeed  profile  closely 
resembles  this  prediction.  The  temperature-driven  profile  displays  a  poor  transmission  in 
the  5  m  channel,  but  a  transmission  nearly  identical  to  the  isospeed  profile  in  the  10  m 
channel,  and  closely  approximating  the  Marsh  &  Schulkin  results  for  the  deeper  channels. 
The  salinity-driven  profile  has  good  transmission  at  5  m  depth,  followed  by  the  poorest 
transmission  quality  at  10  m  depth  and  transmission  similar  to  the  temperature-driven 
profile  in  the  deeper  channels. 


Channel  Depth  (m) 


Figure  10.  Impact  of  channel  depth 


H.  BOTTOM  SEDIMENT  MATERIAL 

The  composition  of  the  ocean  bottom  is  an  important  parameter  since  sound 
energy  interacts  with  the  seafloor.  Ocean  bottom  sediment  materials  are  often  modeled  as 
fluids,  so  they  only  support  compressional  sound  waves.  This  is  an  acceptable 
approximation  since  the  shear  speed  of  the  sediment  in  usually  considerably  less  than  that 
of  a  solid.  High-frequency  sound  energy  does  not  deeply  penetrate  the  seabed,  so  it  is 
acceptable  to  assume  that  the  entire  seafloor  is  homogenous,  neglecting  the  underlying 
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strata  below  the  sediment.15  Parametric  analysis  of  transmission  quality  for  different 
bottom  sediment  materials  indicates  which  geoacoustic  property  is  more  important  in 
determining  the  transmission  absorption  characteristics  of  the  material.  Table  3  shows 
the  geoacoustic  properties  of  the  most  common  bottom  sediment  materials.  ptJ  pw  is  the 
ratio  of  the  material  density  to  seawater  density,  cp  is  the  compressional  wave  speed,  cp/cw 
is  the  ratio  of  material  compressional  wave  speed  to  seawater  compressional  wave  speed, 
and  ap  is  the  compressional  wave  attenuation. 


Bottom  Type 

Porosity 

(%) 

Pi/  Pw 

Cp 

(m/s) 

£// 

ap 

(dB/X„) 

Sand 

45 

1.9 

1650 

1.10 

0.8 

Clay 

70 

1.5 

1500 

1.00 

0.2 

Silt 

55 

1.7 

1575 

1.05 

1 

Gravel 

35 

2.0 

1800 

1.20 

0.6 

Table  3.  Geoacoustic  properties  of  channel  bottom  sediment  materials16 


m 

I 


Bottom  Sediment  Material 


Figure  1 1 .  Impact  of  bottom  sediment  material 


15  Jensen,  38. 

16  Ibid. 
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Figure  11  is  the  summary  of  the  tests  done  for  different  bottom  sediment 
materials.  The  figure  shows  that  the  temperature-driven  and  salinity-driven  environments 
exhibit  similar  trends  regarding  the  order  of  the  strength  of  the  signals  at  the  receiver.  For 
both  cases,  the  silt  sediment  gave  the  strongest  signal,  followed  by  sand,  gravel  and  clay. 
The  only  geoacoustic  parameter  in  which  the  silt  has  the  greatest  magnitude  is  the 
compressional  wave  attenuation,  and  the  order  in  magnitudes  for  this  parameter  reflects 
the  order  of  transmission  strength,  so  the  value  of  the  compressional  wave  attenuation  is 
hypothesized  to  be  the  dominating  factor  in  attenuating  the  signal. 

For  an  isospeed  environment,  the  value  of  the  maximum  path  strength  does  not 
change  since  the  maximum  amplitude  is  the  result  of  a  direct  path  between  the  source  and 
the  receiver,  which  is  unaffected  by  the  characteristics  of  the  bottom  sediment.  As' 
expected,  the  secondary  arrivals  are  observed  to  follow  the  same  trend  as  the  maximum 
arrivals  in  the  salinity  and  temperature-driven  cases. 

I.  TRANSMITTER  AND  RECEIVER  GEOMETRY 

Nine  different  combinations  of  receiver  and  transmitter  height  in  the  channel  were 
tested  for  each  sound-speed  profile.  The  transmitter  and  receiver  were  placed  at  heights 
of  1,  2  and  4  m  above  the  seabed.  Figure  12  shows  the  amplitude  at  the  receiver  for 
transmissions  in  the  standard  channel  with  a  salinity-driven  environment.  The  figure 
shows  that  eight  out  of  the  nine  combinations  display  very  similar  transmission  qualities. 
The  transmission  with  both  transmitter  and  receiver  at  4  m  is  hypothesized  to  have  higher 
amplitude  due  to  the  presence  of  a  strong  single  surface  bounce  path. 
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Transmitter  Depth  (m)  Receiver  Depth  (m) 


Figure  12.  Amplitude  for  different  transmitter  and  receiver  geometries  for  a 

salinity-driven  environment 

Figure  13  shows  the  amplitude  at  the  receiver  for  transmissions  in  the  standard 
channel  with  a  temperature-driven  profile.  As  in  the  salinity-driven  profile,  most  of  the 
transmissions  are  observed  to  have  similar  amplitudes.  In  the  temperature-driven  profile, 
however,  it  is  observed  that  transmissions  with  receiver  and  transmitter  at  the  same  height 
are  stronger  than  those  with  dissimilar  heights.  The  transmission  with  both  transmitter 
and  receiver  at  4  m  is  hypothesized  to  have  the  highest  amplitude  due  to  a  strong  single 
bottom  bounce  path.  This  amplitude  is  slightly  less  than  the  maximum  amplitude  in  the 
salinity-driven  profile  because  of  losses  in  the  bottom  boundary. 
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Transmitter  Depth  (m)  Receiver  Depth  (m) 

Figure  13.  Amplitude  for  different  transmitter  and  receiver  geometries  for  a 

temperature-driven  environment 

Figure  14  shows  the  amplitude  at  the  receiver  for  transmissions  in  the  standard 
channel  with  an  isospeed  profile.  The  figure  shows  no  change  in  the  maximum  amplitude 
with  changes  in  transmitter  and  receiver  geometry.  This  is  because  the  maximum 
amplitude  in  this  case  is  from  an  unrefracted  direct  path  between  source  and  receiver,  and 
no  changes  in  transmitter-receiver  geometry  affect  that  straight-line  path. 
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Transmitter  Depth  (m)  Receiver  Depth  (m) 

Figure  14.  Amplitude  for  different  transmitter  and  receiver  geometries  for  an 

isospeed  environment 
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V.  MARITIME  ENVIRONMENT  CHARACTERIZATION 


A.  ST.  ANDREW  BAY,  PANAMA  CITY,  FLORIDA 


The  first  maritime  environment  considered  is  St.  Andrew  Bay,  Florida,  where 
transmission  experiments  are  taking  place  in  June  2007.  Sound-speed  profile  data  were 
obtained  during  testing  here  in  July  2003.  Conductivity-Temperature-Depth  (CTD) 
measurements  were  performed  through  the  channel  to  determine  the  sound-speed  profile. 
Figure  15  shows  the  bathymetry  at  the  site. 


St.  Andrew  Bay,  Panama  City,  FL 
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Figure  15.  Bathymetric  chart  of  St.  Andrew  Bay,  Panama  City,  FL 
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The  average  channel  depth  for  this  location  is  10  m,  and  the  desired  transmission 
range  is  500  m.  Planned  transmissions  in  this  channel  are  in  the  9-14  kHz  band,  with  the 
model  assuming  a  14-kHz  transmission  since  the  upper  band  displays  the  highest  amount 
of  transmission  loss  due  to  frequency-dependence.  The  2003  test  revealed  that  the  bay  is 
characterized  by  the  presence  of  warmer,  fresher,  less  conductive,  and  less  dense  water 
overlaying  colder,  more  saline,  more  conductive,  denser  water  due  to  fresh  water  from 
rivers  and  rain  run-off  with  variations  in  the  profde  due  to  the  interaction  of  these  two 
layers.  The  sound-speed  profiles  were  compiled  over  a  span  of  7  days.  The  average  SSP 
for  three  different  days  are  illustrated  in  Figure  16. 


Figure  16.  July  2003  average  sound-speed  profiles  in  the  shallow-water 

channel  in  St.  Andrew  Bay,  Panama  City,  FL 

The  figure  shows  that  all  profiles  are  salinity-driven,  with  sound  speed  increasing 
with  depth,  resulting  in  upward-refracting  propagation.  The  test  also  found  that  the 
bottom  material  composition  was  a  mud/silt  combination. 
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Figure  17.  Amplitude  at  the  receiver  for  different  SSP  measured  in  St. 

Andrew  Bay,  FL 


Figure  17  shows  the  maximum  amplitude  for  each  of  the  profiles  shown  in  Figure 
16.  Figure  17  shows  that  there  can  be  drastic  changes  to  the  transmission  characteristics 
of  a  shallow  water  channel  for  slight  changes  in  the  profile  and  for  small  lapses  in  time. 
The  sound-speed  profile  for  25  July  causes  two  families  of  multipaths  to  be  formed, 
which  creates  two  modes  of  arrivals  with  similar  magnitudes  in  very  short  periods  of 
time.  These  multiple  arrivals  have  small  magnitude,  however,  and  do  not  constructively 
interfere  with  one  another.  The  rest  of  the  impulse  response  plot  is  similar  to  those  of  the 
other  days.  The  impulse  response  plots  for  30  July  and  31  July  have  similar  impulse 
distribution,  and  the  higher  amplitude  for  30  July  is  hypothesized  to  be  due  to  the 
convergence  at  the  receiver  of  multiple  rays  experiencing  a  single  surface  reflection. 
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B. 


CHESAPEAKE  BAY,  LITTLE  CREEK,  VIRGINIA 


The  second  maritime  environment  is  Chesapeake  Bay,  Virginia,  where 
transmission  experiments  are  taking  place  in  June  2007.  Figure  18  shows  the  bathymetry 
at  the  site. 
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Figure  18.  Bathymetric  chart  of  Chesapeake  Bay,  Little  Creek,  VA 
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The  typical  channel  depth  for  this  location  is  5  m,  and  the  desired  transmission 
range  is  1000  m.  As  in  St.  Andrew  Bay,  planned  transmissions  in  this  channel  are  in  the 
9-14  kHz  band,  with  the  model  assuming  a  14-kHz  transmission  since  the  higher 
frequencies  are  prone  to  greater  transmission  loss.  The  sound-speed  profiles  for  this 
analysis  were  the  slowest,  fastest,  and  average  SSP  measured  during  testing.  These 
sound-speed  profiles  are  illustrated  in  Figure  19. 
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Figure  19.  Measured  sound-speed  profiles  in  the  shallow-water  channel  in 

Chesapeake  Bay,  Little  Creek,  VA 

The  figure  shows  that  the  slowest  measured  SSP  is  linear  and  very  slightly 
upward  refracting,  while  the  fastest  measured  SSP  is  isospeed  from  the  surface  to  a  2  m 
depth,  then  strongly  downward  refracting.  The  average  SSP  is  slightly  downward 
refracting  from  the  top  to  a  3  m  depth,  then  isospeed  until  the  bottom. 
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Amplitude  at  the  receiver  for  a  1  W  Source  for  SSP  in  Chesapeake  Bay,  VA 
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Figure  20.  Amplitude  at  the  receiver  for  the  measured  sound-speed  profiles 

in  Chesapeake  Bay,  Little  Creek,  VA 


Figure  20  shows  the  maximum  amplitude  for  each  of  the  profiles  shown  in  Figure 
19.  Figure  20  shows  that  the  slowest  measured  SSP  has  the  greatest  amplitude  at  the 
receiver,  and  that  the  fastest  measured  SSP  has  the  smallest  amplitude  at  the  receiver. 
Analysis  shows  that  the  fastest  SSP  concentrates  most  of  the  beams  and  energy  in  the 
isospeed  part  of  the  channel,  with  only  a  fraction  reaching  the  bottom.  This  is 
hypothesized  to  be  the  cause  of  the  small  amplitude  at  the  receiver.  The  slowest  SSP  is 
very  slightly  upward  refracting  and  thus  does  not  severely  affect  the  energy  distribution 
in  the  water  channel.  The  average  SSP  develops  a  duct  at  2.75  m  that  concentrates  most 
of  the  sound  energy  at  the  bottom  of  the  channel.  The  duct  in  the  average  SSP  case 
supports  a  direct  path  from  source  to  receiver  which  curiously  gives  a  lower  amplitude 
than  the  surface-reflected  paths  of  the  slowest-measured  SSP  case.  These  amplitudes 
were  confirmed  with  Bellhop  runs  using  very  dense  ray  fans. 
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VI.  CONCLUSIONS 


A.  FINDINGS 

The  Bellhop  propagation  model  is  a  useful  and  intuitive  tool  for  assessing  the 
propagation  characteristics  of  shallow  water  channels.  Ray-trace  plots  display  the 
distribution  of  sound  energy  in  the  water  channel,  while  eigenray  plots  identify  the  paths 
that  transport  energy  from  the  transmitter  to  a  specific  receiver  location.  The  impulse 
response  plots  show  the  multipath  time  dispersion  of  the  acoustic  signal  arising  from  the 
difference  in  pathlength  traveled  by  each  ray. 

In  the  shallow  waters  of  coastal  regions,  the  velocity  profile  tends  to  be  irregular 
and  unpredictable,  with  little  temporal  or  spatial  stability.  The  shallow  water  profile  is 
complicated  by  the  effects  of  salinity  changes  caused  by  nearby  sources  of  fresh  water. 
Characterizing  real  sound  velocity  profiles  results  in  more  complicated  ray-trace  and 
impulse  response  plots  due  to  the  irregularities  in  the  SSP.  Careful  examination  of  the 
channel  conditions  is  imperative  before  system  deployment  to  ensure  successful 
communications.  The  propagation  model  implemented  here  is  a  useful  tool  for  mission 
planning. 

The  Marsh  &  Schulkin  equation  closely  matches  the  results  obtained  from  the 
Bellhop  algorithm. 

Parametric  analysis  shows  that  transmission  range  is  the  factor  giving  the  largest 
impact  on  signal  strength.  Channel  depth,  transmission  frequency  and  bottom  sediment 
material  also  affect  the  quality  of  transmission  but  their  effect  is  less  than  that  of  range. 
Transmitter  and  receiver  geometry  can  have  a  discernible  effect  on  transmissions  under 
specific  conditions.  A  library  of  cases  is  included  as  Appendix  B  to  illustrate  how 
changes  in  each  parameter  affect  the  ray  structure  and  impulse  response  of  the 
transmission. 
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B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

There  is  a  need  to  perform  experiments  to  validate  and  tune  the  Bellhop  model. 
This  experiments  need  to  occur  in  a  controlled  environment  in  which  each  parameter  can 
be  varied  individually. 

Additionally,  the  assumption  here  that  the  sea  surface  is  a  perfect  reflector  should 
be  examined.  An  additional  parametric  study  involving  sea  states  greater  than  sea-state  0 
should  be  undertaken  to  examine  the  impact  of  the  sea  surface  on  propagation. 

Once  this  work  is  completed,  the  Bellhop  model  may  be  used  to  predict 
performance  of  acoustic  communications  in  maritime  water  channels. 
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APPENDIX  A.  MATLAB  CODES 


A.  COMPOSITE  FIGURE  SUBROUTINE 

%Composite  Figure  Subroutine 

%Open  Sound-Speed  Profile  (SSP)  from  text  file  and  plot  it 

fid=fopen('filename.txf); 

xy=textscan(fid,'%f  %f ); 

x=xy{2}; 

y=xy{l}; 

subplot(3,4,[l  5]) 
plot(x,y  ,'k','Line  W  idth'  ,2) 
set(gca,'YDir','reverse') 
xlabel('Sound  Speed  (m/s)') 
ylabel('Depth  (m)') 
title('Sound  Speed  Profile') 
axis([  1500.0  1550.0  0.0  10.0]) 

%Import  Eigenray  plot 
subplot(3,4,[2  8]) 
hold  on 

plotray('filename') 


39 


%Import  Impulse  Response  plot 
subplot(3,4,[9  12]) 
hold  on 
filenameA 


B.  IMPULSE  RESPONSE  PLOT  SUBROUTINE 


%Impulse  Response  Plot  Subroutine 


%Open  Amplitude-delay  arrival  file 
fid=fopen('filename .  arr') ; 

%Ignore  the  first  six  lines  of  the  amplitude-delay  arrival  file  code  (which  display 
frequency,  receiver  depth,  source  depth,  transmission  range  and  number  of  arrivals) 

for  i=l:6 


tline=fgetl(fid); 

end 

xy=textscan(fid,'%f  %f  %f  %f  %f  %d  %d’); 

x=xy{3}; 

z=xy{l}; 

y=170.8+20*logl0(z); 
stem(x,y,'b.') 
xlabel('Time  (s)') 
ylabel(' Amplitude  in  dB') 

title('Channel  Impulse  Response  for  a  SL  of  1  W') 
axis([0  1  90.0  140.0]) 
axis  'auto  x' 
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APPENDIX  B.  COMPOSITE  FIGURES 


A.  PARAMETRIC  ANALYSIS 


Bellhop  Analysis,  Case  10.500.1. 1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1. 1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.70.Sand.-2 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

-2 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


43 


Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1. 1.70.Sand.-l 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

-1 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


0.332 
Time  (s) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1. 1.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

0  (Isospeed) 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


0i - • - < — — . - . - 

1  - 
2  - 

3  - 

4  - 

5  - 

6  - 

7  - 

8  - 
9  - 

loj - ' - '—I — ' - 1 - 

1500  1510  1520  1530  1540  1550 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.70.Sand.+l 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

+1 

BELLHOP-  Ray  Trace 


Range  (m) 
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Bellhop  Analysis,  Case  10.500.1. 1.70.Sand.+2 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

+2 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Time  (s) 
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Bellhop  Analysis,  Case  10.500.1. 1.70.Clay.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Clay 

SSP  Gradient 

T  emperature-Driven 
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Amplitude  in  dB  ^  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.70.Clay.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Clay 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


0.3305  0.331 

Time  (s) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.70.Clay.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Clay 

SSP  Gradient 

Isospeed 

Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


1500  1510  1520  1530  1540  1550 
Sound  Speed  (m/s) 


250 
Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.70.Silt.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Silt 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 


Range  (m) 


51 


Amplitude  in  dB  Depth  ^  Depth 


Bellhop  Analysis,  Case  10.500.1. 1.70.Silt.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Silt 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


0.3305  0.331 

Time  (s) 


52 


Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.70.Silt.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Silt 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


1500  1510  1520  1530  1540  1550 
Sound  Speed  (m/s) 


250 
Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.70.Gravel.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Gravel 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


1 

! 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W- 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.70.Gravel.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Gravel 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


250 
Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.70.Gravel.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Gravel 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


0i - . - . — — . - . - 

1  - 
2  - 

3  - 

4  - 

5  - 

6  - 

7  - 

8  - 
9  - 

1&X)  1510  1520  1530  1540  1550 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.65.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

65 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Time  (s) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.65.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

65 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1. 1.65.Sand.O 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

65 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.60.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

60 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1. 1.60.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

60 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W' 


250  300 

Range  (m) 
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Amplitude  in  dB  ^  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.60.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

60 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


0, - . - i — — . - « - 

1  - 
2  - 

3  - 

4  - 

5  - 

6  - 

7  - 

8  - 
9  - 

10l - , - >— I — , - , - 

1500  1510  1520  1530  1540  1550 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.55.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

55 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.55.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

55 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 


Range  (m) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1. 1.55.Sand.O 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

55 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


0i - • - < — — . - . - 

1  - 
2  - 

3  - 

4  - 

5  - 

6  - 

7  - 

8  - 
9  - 

loj - ' - '—I — ' - 1 - 

1500  1510  1520  1530  1540  1550 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.50.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

50 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

0 

1 

2 

3 


t: 

7 

8 
9 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.50.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

50 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Channel  Impulse  Response  for  a  SL  of  1  W 


Time  (s) 


67 


Amplitude  in  dB  ^  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.50.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

50 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


0i - • - < — — . - . - 

1  - 
2  - 

3  - 

4  - 

5  - 

6  - 

7  - 

8  - 
9  - 

loj - ' - '—I — ' - 1 - 

1500  1510  1520  1530  1540  1550 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.45.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

45 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.45.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

45 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.45.Sand.O 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

45 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1. 1.40.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

40 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range (m) 
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Amplitude  in  dB  ^  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.40.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

40 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.500.1. 1.40.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

40 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


^OO  1510  1520  1530  1540  1550 
Sound  Speed  (m/s) 


100  150  200  250  300 

Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


0.331  0.332 

Time  (s) 


74 


Bellhop  Analysis,  Case  10.500.1.2.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


1 

! 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W- 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1.2.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1.2.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


1500  1510  1520  1530  1540  1550 
Sound  Speed  (m/s) 


100  150  200  251 

Range 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.1.4.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


1 

! 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W- 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1.4.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.1.4.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Channel  Impulse  Response  for  a  SL  of  1  W' 


BELLHOP-  Ray  Trace 


50  100  150  200  250  300  350  400  450  500 

Range  (m) 


80 


Bellhop  Analysis,  Case  10.500.2.1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.2.1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.2.1.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.2.2.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


1 

! 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W- 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.2.2.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200  251 

Range 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.2.2.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Channel  Impulse  Response  for  a  SL  of  1  W' 


BELLHOP-  Ray  Trace 


50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Bellhop  Analysis,  Case  10.500.2.4.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


1 

! 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W- 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.2.4.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.2.4.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

2 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


0, - 1 - « — — . - . - 

1  - 
2  - 

3  - 

4  - 

5  - 

6  - 

7  - 

8  - 
9  - 

loj - ' - '—I — ' - 1 - 

1500  1510  1520  1530  1540  1550 


2 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.500.4.1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


90 


Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.4.1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W' 


250  300 

Range  (m) 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.4.1.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Channel  Impulse  Response  for  a  SL  of  1  W' 


BELLHOP-  Ray  Trace 


50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Bellhop  Analysis,  Case  10.500.4.2.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.4.2.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


94 


Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.4.2.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

2 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


1500  1510  1520  1530  1540  1550 
Sound  Speed  (m/s) 


100  150  200 

Channel  Impulse  Response  for  a  SL  of  1  W' 


250  300 

Range  (m) 
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Bellhop  Analysis,  Case  10.500.4.4.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


■§  90  - 
|  70- 

-  0.328  0.329 


_ I _ 

0.33 


Time  (s) 


_ I _ 

0.331 


_ I _ 

0.332 


_ I _ 

0.333 


_ I 

0.334 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.4.4.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  ^  Depth  ^ 


Bellhop  Analysis,  Case  10.500.4.4.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

4 

Receiver  Height  (m) 

4 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.250.1. 1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

250 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.250.1. 1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

250 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.250.1. 1.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

250 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


1500  1510  1520  1530  1540  1550 
Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  10.1000.1. 1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

1000 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


-§  90  - 
|  70- 

0.6554 


i  It  l.ff  I 

0.6556 


0.656  0.6562 

Time  (s) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.1000.1. 1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

1000 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

0  100  200  300  400  500  600  700  800  900  1000 

Range (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  10.1000.1. 1.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

1000 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


0  100  200  300  400  500  600  700  800  900  1000 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


1500  1510  1520  1530  1540  1550 
Sound  Speed  (m/s) 


100  200  300  400  500 

Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Bellhop  Analysis,  Case  5.500.1. 1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

5 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


-§  90 

I  70 


0*3275 


0.3278 
Time  (s) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  5.500.1. 1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

5 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  5.500.1. 1.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

5 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 
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Amplitude  in  dB  ^  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  20.500.1. 1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

20 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  20.500.1. 1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

20 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 


0.34 
Time  (s) 
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Amplitude  in  dB  ^  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  20.500.1. 1.70.Sand.0 


Parameter 

Magnitude 

Channel  Depth  (m) 

20 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range (m) 


90  - 
70  - 

0325  0.33 


I 

0.335 


_ I _ 

0.34 


Time  (s) 


Ji 


_ I _ 

0.35 


_ I 

0.355 
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Amplitude  in  dB  ^  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  40.500.1. 1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

40 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

T  emperature-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  40.500.1. 1.70.Sand.Sal 


Parameter 

Magnitude 

Channel  Depth  (m) 

40 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Salinity-Driven 

BELLHOP-  Ray  Trace 


Range  (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  40.500.1. 1.70.Sand.Temp 


Parameter 

Magnitude 

Channel  Depth  (m) 

40 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

70 

Bottom  Sediment  Material 

Sand 

SSP  Gradient 

Isospeed 

BELLHOP-  Ray  Trace 


Range  (m) 


113 


Amplitude  in  dB  Depth  (m)  Depth  (m) 


B. 


ST  ANDREW  BAY,  PANAMA  CITY,  FL 


Bellhop  Analysis,  Case  SAB.10.500.1.1.14.Silt.0725 


Parameter 

Magnitude 

Location 

St.  Andrew  Bay 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

14 

Bottom  Sediment  Material 

Silt 

SSP  Gradient 

25  July  2003 

BELLHOP-  Ray  Trace  Plot 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


0.325  0.326  0.327  0.328  0.329  0.33  0.331  0.332  0.333  0.334  0.335 

Time  (s) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  SAB.10.500.1.1.14.Silt.0730 


Parameter 

Magnitude 

Location 

St.  Andrew  Bay 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

500 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

14 

Bottom  Sediment  Material 

Silt 

SSP  Gradient 

30  July  2003 

BELLHOP-  Ray  Trace  Plot 


Range  (m) 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


Bellhop  Analysis,  Case  SAB.10.1000.1.1.14.Silt.0731 


Parameter 

Magnitude 

Location 

St.  Andrew  Bay 

Channel  Depth  (m) 

10 

Transmission  Range  (m) 

1000 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

14 

Bottom  Sediment  Material 

Silt 

SSP  Gradient 

31  July  2003 

BELLHOP-  Ray  Trace  Plot 


0  50  100  150  200  250  300  350  400  450  500 

Range  (m) 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  (m)  Depth  (m) 


C. 


CHESAPEAKE  BAY,  LITTLE  CREEK,  VA 
Bellhop  Analysis,  Case  CB.5.1000.1.1.14.Sand.SM 


Parameter 

Magnitude 

Location 

Chesapeake  Bay 

Channel  Depth  (m) 

5 

Transmission  Range  (m) 

1000 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

14 

Bottom  Sediment  Material 

Sand 

SSP 

Slowest  Measured 

BELLHOP-  Ray  Trace  Plot 


Range  (m) 


Sound  Speed  Profile  BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s)  Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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Amplitude  in  dB  Depth  (m) 


Bellhop  Analysis,  Case  CB.5.1000.1.1.14.Sand.FM 


Parameter 

Magnitude 

Location 

Chesapeake  Bay 

Channel  Depth  (m) 

5 

Transmission  Range  (m) 

1000 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

14 

Bottom  Sediment  Material 

Sand 

SSP 

Fastest  Measured 

BELLHOP-  Ray  Trace  Plot 


400  500  600 

Range  (m) 


90  - 
70  - 

0*6528  0.653 


0.6532  0.6534 


I  T 

0.6536  0.6538 

Time  (s) 


Jt _ L 

0.654 


JU _ l_t _ I _ I 

0.6542  0.6544  0.6546  0.6548 
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Bellhop  Analysis,  Case  CB.5.1000.1.1.14.Sand.AVG 


Parameter 

Magnitude 

Location 

Chesapeake  Bay 

Channel  Depth  (m) 

5 

Transmission  Range  (m) 

1000 

Transmitter  Height  (m) 

1 

Receiver  Height  (m) 

1 

Frequency  (kHz) 

14 

Bottom  Sediment  Material 

Sand 

SSP 

Average 

BELLHOP-  Ray  Trace  Plot 


Range  (m) 


1 

s 


Sound  Speed  Profile 


BELLHOP-  Eigenray  Plot 


Sound  Speed  (m/s) 


Range  (m) 

Channel  Impulse  Response  for  a  SL  of  1  W 
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APPENDIX  C.  SUMMARY  TABLES 


A.  SSP 


SSP 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

T  emperature-Driven 

4.38xl0'4 

103.6 

-1.1 

Salinity-Driven 

2.13xl0'4 

97.4 

-7.1 

Linear  with  a  -2  Slope 

1.91xl0'3 

116.4 

11.1 

Linear  with  a  -1  Slope 

2.10xlCT4 

97.2 

-7.2 

Isospeed 

5. 04x1  O'4 

104.8 

— 

Linear  with  a  +1  Slope 

1.63x1  O'4 

95.0 

-9.3 

Linear  with  a  +2  Slope 

5.21xl0'5 

85.1 

-18.8 

Table  4.  Amplitude  response  for  different  sound-speed  profiles  (SSP) 


B.  TRANSMITTER  FREQUENCY 

1.  Salinity-Driven  Environment 


Frequency 

(kHz) 

1 P(r,  z)\!  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

40 

3. 29x1 0'4 

101.1 

3.8 

45 

2.97x1 0'4 

100.3 

2.9 

50 

3.23xl0~4 

101.0 

3.7 

55 

2.88xl0'4 

100.0 

2.7 

60 

2.59xl0'4 

99.1 

1.7 

65 

2.34x1 0'4 

98.2 

0.8 

70 

2.13xl0'4 

97.4 

— 

Table  5.  Amplitude  response  in  a  salinity-driven  environment  for  different  transmission 

frequencies 
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2.  Temperature-Driven  Environment 


Frequency 

(kHz) 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

40 

8.51xl0'4 

109.4 

5.6 

45 

7.51xl0'4 

108.3 

4.5 

50 

6.65x1  O'4 

107.3 

3.5 

55 

5. 92x1  O'4 

106.2 

2.6 

60 

5. 32x1  O'4 

105.3 

1.7 

65 

4.81xl0'4 

104.4 

0.8 

70 

4.3  8x1  O'4 

103.6 

— 

Table  6.  Amplitude  response  in  a  temperature-driven  environment  for  different 

transmission  frequencies 


3.  Isospeed  Environment 


Frequency 

(kHz) 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

40 

9.80xl0~4 

110.6 

5.6 

45 

8.64x1 0'4 

109.5 

4.5 

50 

7.66xl0'4 

108.5 

3.5 

55 

6.82xl0'4 

107.5 

2.6 

60 

6.12xl0'4 

106.5 

1.7 

65 

5.54xl0'4 

105.7 

0.8 

70 

5. 04x1 0'4 

104.8 

— 

Table  7.  Amplitude  response  in  an  isospeed  environment  for  different  transmission 

frequencies 
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C.  TRANSMISSION  RANGE 

1.  Salinity-Driven  Environment 


Transmission  Range 
(m) 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

250 

2.50xl0'3 

118.8 

21.9 

500 

2.13xl0'4 

97.4 

— 

1000 

1.03xl0'5 

71.1 

-27.0 

Table  8.  Amplitude  response  in  a  salinity-driven  environment  for  different  transmission 

ranges 


2.  Temperature-Driven  Environment 


Transmission  Range 
(m) 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

250 

5. 46x1  O’3 

125.5 

21.2 

500 

4.38xl0'4 

103.6 

— 

1000 

2.06xl0'5 

11 A 

-25.6 

Table  9.  Amplitude  response  in  a  temperature-driven  environment  for  different 

transmission  ranges 


3.  Isospeed  Environment 


Transmission  Range 
(m) 

\p(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

250 

2.01xl0'3 

116.9 

11.5 

500 

5. 04x1 0'4 

104.8 

— 

1000 

6.36xl0'5 

86.9 

-17.1 

Table  10.  Amplitude  response  in  an  isospeed  environment  for  different  transmission  ranges 
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D. 


CHANNEL  DEPTH 


1.  Salinity-Driven  Environment 


Channel  Depth 

(m) 

1 P(r,  z)  |/  [ po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

5 

7.09x1  O'4 

107.8 

10.7 

10 

2.13xl0'4 

97.4 

— 

20 

6.22x1  O'4 

106.7 

9.5 

40 

5.09xl0-4 

104.9 

7.7 

Table  11.  Amplitude  response  in  a  salinity-driven  environment  with  different  channel 

depths 


2.  Temperature-Driven  Environment 


Channel  Depth 

(m) 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

5 

1.07xl0'5 

71.4 

-31.1 

10 

4.3  8x1 0'4 

103.6 

— 

20 

6.52x1 0'4 

107.1 

3.4 

40 

5.04x1 0'4 

104.8 

1.2 

Table  12.  Amplitude  response  in  a  temperature-driven  environment  with  different  channel 

depths 


3.  Isospeed  Environment 


Channel  Depth 

(m) 

1 P(r,  z)\l  \p0\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

5 

5. 04x1 0'4 

104.8 

0.0 

10 

5. 04x1 0'4 

104.8 

— 

20 

5.04x1 0'4 

104.8 

0.0 

40 

5. 04x1 0'4 

104.8 

0.0 

Table  13.  Amplitude  response  in  an  isospeed  environment  with  different  channel  depths 
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E. 


BOTTOM  SEDIMENT  MATERIAL 


1.  Salinity-Driven  Environment 


Bottom  Type 

1 Pfr,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

Sand 

2.13xl0'4 

97.4 

— 

Clay 

1.37xl0'4 

93.5 

-4.0 

Silt 

1.78x1  O'4 

95.8 

-1.6 

Gravel 

2.29x1  O'4 

98.0 

0.6 

Table  14.  Amplitude  response  in  a  salinity-driven  environment  with  different  bottom 

sediment  materials 


2.  Temperature-Driven  Environment 


Bottom  Type 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

Sand 

4.38xl0'4 

103.6 

— 

Clay 

1.82x1  O'4 

96.0 

-7.3 

Silt 

3. 06x1  O'4 

100.5 

-3.0 

Gravel 

5.05xl0"4 

104.9 

1.2 

Table  15.  Amplitude  response  in  a  temperature-driven  environment  with  different  bottom 

sediment  materials 


3.  Isospeed  Environment 


Bottom  Type 

1 P(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

Sand 

5.04x1  O’4 

104.8 

— 

Clay 

5.04x1 0'4 

104.8 

0.0 

Silt 

5. 04x1 0'4 

104.8 

0.0 

Gravel 

5.04x1 0'4 

104.8 

0.0 

Table  16.  Amplitude  response  in  an  isospeed  environment  with  different  bottom  sediment 

materials 


125 


F. 


TRANSMITTER  AND  RECEIVER  HEIGHT 


1.  Salinity-Driven  Environment 


Transmitter 
Height  (m) 

Receiver 
Height  (m) 

\p(r,  z)\I  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

1 

1 

2.13xl0'4 

97.4 

— 

1 

2 

2.26x1 0"4 

97.9 

0.5 

1 

4 

1.88xl0'4 

96.3 

-1.1 

2 

1 

2.05xl0'5 

97.0 

-0.4 

2 

2 

1.89x1  O'4 

96.3 

-1.1 

2 

4 

1.90x1  O'4 

96.4 

-1.1 

4 

1 

1.87xl0'4 

96.2 

-1.2 

4 

2 

1.85xl0'4 

96.1 

-1.3 

4 

4 

1.56xl0'3 

114.7 

17.7 

Table  17.  Amplitude  response  for  different  transmitter  and  receiver  depth  combinations  in  a 

salinity-driven  environment 


2.  Temperature-Driven  Environment 


Transmitter 
Height  (m) 

Receiver 
Height  (m) 

\p(r,  z)\l  \po\ 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

1 

1 

4.38xl0"4 

103.6 

— 

1 

2 

2.50xl0'4 

98.8 

-4.7 

1 

4 

2.43x1 0"4 

98.5 

-4.9 

2 

1 

2. 60x1  O'4 

99.1 

-4.3 

2 

2 

3.30xl0'4 

101.2 

-2.3 

2 

4 

2. 10x1  O'4 

97.2 

-6.1 

4 

1 

2.50xl0'4 

98.8 

-4.7 

4 

2 

2.12xl0'3 

97.3 

-6.1 

4 

4 

1.42x1 0'3 

113.8 

9.9 

Table  18.  Amplitude  response  for  different  transmitter  and  receiver  depth  combinations  in  a 

temperature-driven  environment 
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3.  Isospeed  Environment 


Transmitter 
Height  (m) 

Receiver 
Height  (m) 

\p(r,  z)  |/  M 

Amplitude 

(dB) 

Baseline  Environment  Deviation 

(%) 

1 

1 

5. 04x1  O'4 

104.8 

— 

1 

2 

5. 04x1  O'4 

104.8 

0.0 

1 

4 

5. 04x1  O'4 

104.8 

0.0 

2 

1 

5. 04x1  O'4 

104.8 

0.0 

2 

2 

5. 04x1  O'4 

104.8 

0.0 

2 

4 

5. 04x1  O'4 

104.8 

0.0 

4 

1 

5. 04x1  O'4 

104.8 

0.0 

4 

2 

5. 04x1  O'4 

104.8 

0.0 

4 

4 

5. 04x1  O'4 

104.8 

0.0 

Table  19.  Amplitude  response  for  different  transmitter  and  receiver  depth  combinations  in 

an  isospeed  environment 
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